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We extend a recent analysis of gravitational perturbations on Dirac-Goto-Nambu type strings, 
membranes and higher dimensional branes. In an arbitrary gauge, it is shown that the relevant first 
order equations governing the displacement vector of the worldsheet and metric perturbation are 
obtainable from a variational principle whose Lagrangian is constructed as a second order perturba- 
tion of the standard Dirac-Goto-Nambu action density. A symplectic current functional is obtained 
as a by-product that is potentially useful for the derivation of conservation laws in particular cir- 
Q\ ' cumstances. 

' I. INTRODUCTION 

o ' 

The purpose of this work is to extend our recent analysis [jjj of the dynamics of small perturbations of the worldsheet 
for a simple, internally structureless Dirac-Goto-Nambu type string, membrane, or higher dimensional brane in a 
pseudo-Riemannian background of arbitrary dimension. The main motivation for this investigation is to prepare the 
way for a study radiation backreaction that has already been undertaken |^,|3| in the technically simpler, but physically 
more speculative context of scalar axion fields. In view of this ultimate objective, our analysis jjj specifically allowed 
for the effect of perturbing gravitational waves, unlike earlier work that was restricted to effectively free perturbations 
of the worldsheet in a fixed flat Q or curved space-time background. The main result of this analysis to date 
has been the demonstration |8) that the divergent part of the gravitational self interaction exactly cancels for a 
Nambu-Goto string in ordinary four-dimensional space-time, contrary to previous claims in the literature ]s|Jlfj[|, but 
qq ' consistently with concurrent work [jll| . 

Most simple treatments of the dynamics of branes use a formalism originally developed by Eisenhart Q , using the 
explicit details of a specifically chosen worldsheet reference scheme at each step. However, particularly when there 
are other physically independent helds (gravitational or otherwise) to be taken into account, it is efficient to use a 
more economical treatment in which the amount of auxiliary mathematical structure involved is reduced. The 
approach we have used |jj],[3) reduces the mathematical paraphernalia to its strict minimum by working exclusively in 
i <~| terms of fields whose tensorial components are defined directly in terms of the background space time coordinates, 
that is, without any reference to special frames or internal coordinates. While particular coordinate systems are 
• genuinely useful for specific applications, the excessive number of indices can lead to the mistakes which have already 

■ been discussed above. 

Although it did not allow for gravitational wave perturbations, the first of the works discussed above |j| went 
further than our more recent generalization Q|. In both cases, the first-order perturbation equations governing the 
dynamical evolution of the relevant fields (a vectorial surface field £ M in the case of ref. || and the surface field 
plus the gravitational wave perturbation for ref. [jjj) were deduced. But for the simpler case of just the surface 
field, it was shown that these dynamical equations could be derived via a variational principle, from a Lagrangian of 
quadratic order in £ p , which was constructed as a second order perturbation of the original action — in this case the 
Dirac-Goto-Nambu action. A by-product of this second-order approach is the construction of a bi-linear symplectic 
surface current satisfying a Noether type conservation law. 

In the present work it is shown how to construct an equivalent second-order Lagrangian, and the corresponding 
bi- linear symplectic current in the more general case which includes gravitational wave perturbations. From this 
Lagrangian, we deduce that the same dynamical equations of motion already found in ref. Q and used in ref. ||. 
However, it will be seen that the gravitational waves perturbations act as a source term for the symplectic current, 
and hence in this more general case the relevant Noetherian surface current no longer satisfies a strict conservation 
law. 

We will use a similar notation scheme as before to describe a p-brane with ap + 1 dimensional worldsheet, in 
n space-time dimensions. As emphasized already, the need to refer any internal worldsheet coordinates, a a (a=0, 
p) say, is avoided by basing the analysis on the first fundamental tensor, with components rf ,v defined with 
respect to ordinary background space-time coordinates = X M (er a ) (/i =1, n). This tensor is determined by 
the background space-time metric g^ v as the projection of the contravariant inverse of the induced metric on the 
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worldsheet, that is, in terms of the internal coordinates (using d a for partial differentiation), as r]^ = r y ah d a X p di > X 1 ' , 
where the induced metric is given by 7 a fc = g lllV d a X p 'd\ > X v '. Contracting any vector with this fundamental tensor has 
the effect of projecting it onto the part tangential to the worldsheet, while the normal part is obtainable by contraction 
with the orthogonal projection tensor H 1V = g^ v — rf iV . We denote the operator of tangiantially projected covariant 
differentiation — the only kind that is well defined for a tensor field with support confined to the worldsheet — by 

V^ 1 = T]^V V , (1) 

where V ' v is the usual operator of covariant differentiation as defined with respect to the Riemannian background 
connection, namely T x aT = g Xp (g p ( a ,T) ~ tS<tt,p) using round brackets for index symmetrisation. The information 
characterizing the various (intrinsic and extrinsic) kinds of curvature associated with the worldsheet embedding is 
provided Jl2| by the second fundamental tensor, as defined in terms the first one by 

K^ p = T^V^r ■ (2) 

Since it is automatically surface tangential, as well as symmetric, with respect to its first two indices but surface 
orthogonal with respect to its last index, contraction with the former provides its only non-identically vanishing trace, 
namely the curvature vector 

KP = ^ V K^. (3) 

The non-trivial dynamical requirement that this vector K p should actually vanish is the tensorially covariant 
expression of the equations of motion for the kind of strings and higher branes that we shall be considering, namely 
those characterised by the Dirac-Goto-Nambu action principle. In this case, the relevant action I is just proportional 
to the geometrically induced measure of the worldsheet, meaning that it is obtained by taking a fixed value for the 
Lagrangian density C in the general expression for a worldsheet action, as given by 

1 = J CdS , (4) 
where dS = \\j\\ 1 / 2 d p+1 a . in terms of internal coordinates. 



II. RELATION BETWEEN LAGRANGIAN AND EULERIAN VARIATIONS 



In order to proceed with the second-order analysis of the action principle, we need to define the relevant pertur- 
bations and in particular the difference between Eulerian and Lagrangian variations. In ref. jjj we only included 
the very much simpler first-order perturbations, while in ref. || second-order perturbations were included in a fixed 
background, that is, no Eulerian variation. In this section, we generalize these second-order variations to include 
Eulerian variations of the metric. The starting point for this is the finite action variation, which takes the form 

Al = SI+ 1 -S 2 l+0{5 3 }, (5) 

where S is an infinitesimal variation operator normalised with respect to a variation parameter with magnitude 
proportional to that of the displacement vector field. 
The Eulerian variation of the metric is given by 

9\iv => g^u = g^v + h^v , (6) 

E 

defined with respect to some predetermined position identification scheme, for example, harmonic coordinates. It is 
possible in the analysis of a smoothly distributed fluid to work entirely with Eulerian, that is, 'fixed point', variations. 
However, for the cases under consideration here (point particles, strings, and higher dimensional branes) where the 
relevant fields are confined to a lower dimensional support manifold that may be displaced by the perturbation, it is 
evident that the concept of an Eulerian variation will generically fail to be well defined, so that it becomes necessary 
to use an approach based on Lagrangian, that is, 'comoving', variations. A similar approach is often also used in 
dealing with smooth fluid media. 

As a prerequisite for this, the relevant Eulerian background variation (^) must first be translated into Lagrangian 
form. It is well known that the Lagrangian variation 
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(7) 



representing the difference between the comoving perturbed metric g^ v and the unperturbed value <7 M „ will deviate 

L 

from its Eulerian analogue by a term that to first order in the displacement £ M will be given just by the corresponding 
Lie derivative of the unperturbed field. In order to go to second (and higher) order it is necessary to take care in 
specifying just how the displacement corresponding to £ M is defined. The natural way to do this [[[2| (as suggested, for 
example, by the Bunting identity jl3| , and the analysis of Boisseau and Letelier jlj ) is to specify the displacement 
^ in terms of the corresponding infinitesimal geodesic. On this basis, according to (12) of ref. @, one finds that JTJ] 
that the corresponding Lagrangian value g^ of the perturbed metric will be given in terms of its Eulerian analogue 

L 

g M „ and the associated Riemannian differentiation operator Vp and curvature TZ^pa by 

g?* - a„u = 2V ( ^) + (v^ p ) Vutp - + o{6 3 } . (8) 

L E E EE E 

Subject to the understanding that the metric perturbation is at most of the same infinitesimal order as the displace- 
ment, h^y = 0{<5}, it then follows that the required Lagrangian variation of the metric will be expressible in terms of 
the unperturbed Riemannian differentiation operator and curvature TZp Vp(J by the ubiquitously applicable formula 

- V = 2V (M 6) + £ p V P V + S^V^f + (V„£ p ) V„£„ - + 0{S 3 } . (9) 

The linear order part given by the first term on the right in this formula is the well known Lie derivative term. The 
novelty here is inclusion of the quadratic order adjustment, whose potential importance is to be emphasised: it is 
likely to be useful in more general physical contexts, not just for treating the particular case of the Dirac-Goto-Nambu 
strings and higher branes that are considered here. 

It is to be remarked that there is a lot of arbitrary gauge freedom in the way one chooses the displacement generating 
field £ M at positions off the worldsheet. However, although derivatives of the displacement field in directions orthogonal 
to the worldsheet may play a role in the intermediate steps of the calculations, they must cancel out (a condition 
which provides a useful check on the algebra) in the physical equations governing the motion of the worldsheet, 
for which only the values of £ M at positions actually on the worldsheet are relevant. In the particular case to be 
considered here, Dirac-Nambu-Goto branes, there is another kind of gauge freedom involved, related to the choice of 
the tangential component of the displacement fields. This is because moving along the worldsheet has no effect on its 
locus, which — in the absence of internal currents and fields — is the only physical structure that matters. Hence, 
in the Dirac-Goto-Nambu case there is no loss of physical generality if, for the purpose of obtaining well behaved 
hyperbolic dynamical equations, one chooses to eliminate this freedom by fixing the gauge on the worldsheet through 
imposing the condition of orthogonality, T)^^ — 0. In the present article, we shall nevertheless impose no restriction 
on the gauge to be used, leaving open the possibility of including an arbitrary tangential component in the choice of 



III. SECOND ORDER VARIATION OF THE ACTION. 

We are now in a position to vary the action (^) in the form (^|). It will be convenient to use an abbreviation scheme 
in which an acute accent indicates differentiation with respect to the relevant variation parameter, so that in the 
particular case of the action we are able to simply write 6T = 1 and 8 2 X = T. We should note that, since the action 
is global, the distinction between Lagrangian and Eulerian variations does not arise at this stage. 

It is evident from (^) that in this scheme the first order Lagrangian variation of the metric will be given by an 
expression of the familiar form 

Sg^ = V + 2V ( ^) (10) 
while the corresponding second order Lagrangian variation of the metric will be given by 

fg^ = V + 2v ( /,) + 2f>v p v + ±h p{l y v) ( p + 2(v/") v„£ - 2(p( a n lipva . (11) 

This last expression can be slightly simplified by ommission of the first term on the right if the perturbing field h^ v 
is prescribed to scale directly with the variation parameter, so that one simply has = and = 0. However 
one can not analogously exclude the presence of higher order contributions in the expansion £ M = £ M + + ... for 
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the a priori unknown, displacement field. At this stage we do not make any simplification, but subsequent sections 
we shall. 

To allow for the effects of the perturbation on the surface measure, whose first order Lagrangian variation will be 
given [jjj in terms of the fundamental tensor rf 11 by an expression of the familiar form 

S(dS)=^rr(Sg^)dE, (12) 

it is convenient to express the first and second order variations of any action integral of the form (|4|) as 

OZ)dS, Z= ( (0 2 £)dT. (13) 



The measure weighted 'diamond-differential' operator <0 |H| is surface generalisation of the 'boldface-differential' 
operator already used by Friedman and Schutz [[16| ) . When applied to the Lagrangian density it gives 

$£ = S£+ 1 -£rr6g^, (14) 

so that at second-order one obtains 

f£ = fZ+ (6Z)rrSg M u + \c{t v r, pa ~ 2r,^)(S gtlv )Sg p(7 + iZgW • (15) 

Provided any internal fields that may be present satisfy the corresponding dynamical equations, which implies that 
they give no contribution to the first-order variation, and provided the only relevant external field is that of the 
background metric g^, then the variation in (|l3| ) will be expressible just in terms of the surface energy- momentum 
density tensor in the standard form (jj 

SX=iJ (Sg^) dS, T^ = 2-^ + Zrf v . (16) 

We shall now restrict our attention to the simple case of a Dirac-Goto-Nambu membrane or string, for which 
Lagrangian density in (^) is just a constant, 

£=-m p+1 , (17) 

where m is a fixed parameter which in natural units will have the dimensions of mass. For topological defect system, 
such as a cosmic string or domain wall, this Kibble mass parameter will have the same order of magnitude as the Higgs 
mass associated with the underlying spontaneous symmetry breaking. In this simple case, the Lagrangian variation 
of the Lagrangian density will be given trivially by SL = and the surface energy-momentum density tensor ( |l6| ) will 
simply be given in terms of the fundamental tensor by the proportionality relation 

= -m p+ V v . (18) 



Under these circumstances it can immediately be seen from the preceeding formulae (|10|) and (|1J) that the integrand 
of the first order variation of the action will be given by 



m (p+D <>£ - \rf v ^u = \rf v K v + . (19) 



The content of the variation principle is that the surface integral of this quantity must be made independent of the 
infinitesimal displacement vector field In other words the dynamical equations are specified as the condition for 
all the displacement dependent terms on the the right hand side of ([l9]) to be contained within a surface current 
divergence, that is, a current which is tangential to the worldsheet. The relevant term would already have 

the form of a surface current divergence as it stands, if the displacement £ M satisfied the worldsheet tangentiality 
condition £ M = rj^ v . More generally, in the non trivial case for which it has an orthogonal part, 11 L V £ M ^ 0, the 
formula V„ L" = —K^ can be used to expand the term in question as 

%&=%(rf u C) -K^, (20) 
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where is the curvature vector, as defined by (||). Since the first term in ( pCf ) has the form of a surface current 
divergence, and as such is irrelevant from the point of view of the variation principle, we can replace ([l9]) by the 
expression 



m 



-(p +1 ) <> C = V" V - K ^ > (21) 



using the equivalence symbol = to indicate equality modulo divergences that are removable by Green's theorem for a 
variation with compact support. At this first order level the metric variation h^ v cannot couple to the displacement. 
It thus becomes evident that the content of the ensuing Dirac Goto Nambu type field equation is simply the well 



known condition that this vector should vanish 



We are now ready to move onto new ground by using (11) and ( |l5| ) to obtain the second-order analogue of ( fl9| ) 
which works out as 



In the manner already demonstrated in the special case j^] for which the gravitational perturbation h^ v is absent, 
the second variational integrand (|2^) will be employable as the Lagrangian density of a variational principle for the 
perturbed dynamical equations at first-order. For this purpose, the only relevant terms are those having homoge- 
neously quadratic dependence on first-order variations. The remaining terms, that is to say the last two terms which 
are homogeneously linear in second-order variations, are of the same form as the linear contribution ( |l9| ) in the first- 
order contribution, which means that as consequence of the zero order field equations their contribution to the action 
integral will just be a constant. Since the present analysis is based on the treatment of the gravitational perturbation 
as a given background field, the the first term in (^2|) will also just be a constant as far as the variational principle is 
concerned 



IV. DERIVATION OF THE PERTURBED DYNAMICAL EQUATIONS 

Therefore, we have argued that dropping the constant contributions (from the first term in ( p2| ) and from the 
second-order variations, which are irrelevant for the application of the variation principle, one can obtain a second- 
order action Z{ 2 } say, (differing from T in ( |l3| ) only by the omitted constant) of the form 



1{2} = / C {2} dS, (23) 
where the relevant second-order Lagrangian is given by 

+ 14 (2h;v„e + vr^h^) + v pa v„ (vf) . (24) 

To apply the variation principle to 1{2} we need to consider a process in which the linearized field variables 
themselves undergo independent infinitesimal 'virtual' variations £ M — > £ M + <5£ M and h^ v — > + 8h^ v . Writing 
£ M = (5£ M and — Sh^, that is, using a grave accent to indicate differentiation with respect to the 'virtual' 
variation parameter (as distinct from differentiation with respect to the 'real' or 'physical' variations for which the 
acute accent is used), then the ensuing action differential will have the form 



Sl {2} = 2 I l {1<1} dS, (25) 

where ,i} is a symmetric bi-linear functional of the independent 'real' (acute accented) and 'virtual' (grave accented) 
variations, that is given by 



rn 



- (p+1) c {hl} = 2(v [P i p ) v„]r+ k (v^) - ^ v n wva ( p i° 



k ( h P v Vut + v ) + \ k i 1 " (r Vp v + e v VpO + W a v - (vr + vr) ■ (26) 
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At the expense of sacrificing the manifest symmetry between 'real' and 'virtual' contributions, the next step is to 
introduce linear functionals Q = Q{£,} and F — F{h] depending respectively on the displacement vector £ M and the 
gravitational perturbation /i M „, so as rewrite this in the more directly applicable form 

- m-(P +1 )£ {M} = £»(Fp - Qfj) + e% + %> , (27) 

where we have used the obvious notation 

Qn = QA&> F p = F^{h}, F fl =F fl {h}. (28) 

The surface current vector J M is a symmetric bi-linear functional of the two variations, Q a is the dynamical functional 
for the surface perturbation, and F a is the equivalent functional for the gravitational perturbation. These expressions 
are given by, 

> = k + ^(ev.r - e%e)^u ± p « (r v + H?) + + r v) > (29) 



F a {h) = lf a ( V"V P V - V„K) + (Kf a - l? c K v - \tf v KM ■ (31 



One can see that the last functional is automatically worldsheet orthogonal, that is, for an arbitrary perturbing field 

rf p FAh} = 0, (32) 

is an identity. 

The requirement of the variation principle is that the surface integral (p5| ) of ( p7|) should be independent of the 
virtual displacement £ M for a variation with compact support. Since the term involving the current J M is a pure 
surface divergence, this requirement is evidently equivalent just to the condition that for a given 'real' gravitational 
perturbation field h pl , the corresponding 'real displacement vector ^ should satisfy a field equation of the form 

Qa = F a . (33) 

The validity of this field equation can be immediately confirmed by checking that when the unperturbed, zero order 
dynamical equations K p = are satisfied, ( |33"| ) does indeed agree with the first-order perturbed dynamical equations 
given by equation (21) of our proceeding work In particular, it can be verified that the backreaction terms F\ a 
and F20 that were defined there by equation (45) add up to a total that matches the gravitational forcing term on 
the right of (p3|), that is, the linear functional defined here by ( |3l| ) is expressible in terms of our previous notation 
as F a — Fiu + F%u- According to (41) of ref. pi, the precise physical interpretation of this is that the quantity 
j~ u = —m^ ?Jrl ^Fu represents the gravitational force density exerted on the brane by a perturbing field h pu . 

V. THE SYMPLECTIC CURRENT 

Having confirmed that the homogeneous quadratic functional 1{2} given by ( p2| ) does indeed act as a Lagrangian for 
the linearly perturbed field equations, we can exploit it for the purpose of derivation of various kinds of conservation 
law. In case where the gravitational perturbation h pv it was shown that the current treatment yielded a Noether 
type identity. This expresses the fact that, although it is not manifest in the (p7|), this bilinear variation term is 
by construction symmetric in the sense of being invariant under interchange of the two independent, mathematically 
equivalent variations involved. This symmetry property is equivalently expressible as the vanishing of the antisym- 
metric quantity obtained by taking the difference between (|27j) and its analogue obtained by swapping the acute and 
grave accents. This identity which reduces, by the symmetry property of partial differentiation, to the form 

= eQ p - CQ P , (34) 

in terms of a symplectic — meaning antisymmetric bilinear — surface current C M = C^jf,^} given by C = — C = 
jn _ jfj, w h ere jm i s obtained from J M by the interchange of 'real' and 'virtual' variations. This gives 

'& = (iT L pa + 2tf pV v a] ) (i p V„£* - & V„r ) , (35) 
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in which it is to be noticed that which the external gravitational perturbation field has cancelled out. Hence, the 
geometrically defined current functional C M thus retains exactly the same form, apart from the omission of the physical 
weighting factor m p+1 in its definition, as the corresponding dimensionally weighted symplectic current C M that was 
obtained in ref. @. The current (ff) can be rewritten in the form 

d" =± pa (e p v^r - e v"r) + 2^ v CT a ^ ) . (36) 

from which it can be seen that, for the case of worldsheet conserving displacements, that is, when £ M and £ M are both 
tangential to the brane, the first term will simply vanish while the second term will have the form of an exact surface 
divergence. This implies that that it will be trivially conserved as an identity, independently of any field equations 
that may be satisfied |l7j . 

Although it retains the same form, the symplectic current now no longer obeys the simple surface current conser- 
vation law that was obtained Q in the absence of gravitational perturbations. £ M is postulated to satisfy a dynamical 
equation of the same form as the first order perturbation equation (^) that must be satisfied by allowing the 
dependence on the displacements to be eliminated from the right hand side of (|34|), and hence one is left with a 
relation of the form 

V^c" = - &h > (37) 

in which, generically, there remains a source term on the right in cases where the forcing effect of gravitational wave 
perturbations is present. 

By the orthogonality property (^) , the trivial tangential parts of £ M and £ M make no contribution to the right hand 
side of ([37|). To see how they affect the left hand side, it is useful to rewrite the formula ( |36| ) for the current in the 
form 

(> =l po (i"W(±i h - h) + ^ v Ka( [v e ] + z^^U 19 ^) , (38) 

in which the last term has the form of an exact surface divergence, so that independently of the field equations it will 
be trivially conserved as an identity. Provided the unperturbed field equations, = are satisfied, the second term 

simply drops out, so the 'reduced' symplectic current vector consisting just of the first term, namely 

s" =L pa (e"V"(i5 w h - h) , (39) 

will satisfy an equation of the same form as the original symplectic current, that is, 

= i»F„ - &h , (40) 

Unlike the original current C M , this 'reduced' current S 11 has the advantage of being gauge independent in the sense 
of being entirely independent of the tangentially projected parts of the vectors £ M and £ M on which it depends. 

In the kind of practical application for which such a Noetherian identity is typically used, namely for exploiting 
underlying symmetries of various kinds, it is usually sufficient to consider only virtual displacements satisfying the 
unforced dynamical equations, which means taking h^v — 0, and thereby setting i*) t = in the right hand side of 
(|37j). A virtual displacement field satisfying this condition will be given automatically by setting £ M — fc M where fc M 
is any solution of the background Killing equations V( M fc„) = 0, since the action of such a translation on any given 
unperturbed solution evidently translates it onto a nearby solution that is geometrically identical and therefore also 
a solution. What this means - for any real physical solution £^ of the displacement perturbation equation (|3^) - is 
that associated with the geometric symmetry generated by fc^ there will be a corresponding generalised momentum 
current 

s»{k} =t pct (fv"(i; n - e*v"cr v n) , (41) 

to obey a divergence equation of the form 

Vj{kr - • (42) 

In view of the orthogonality property (|3^) the source term on the right will evidently vanish in the case for which 
fc M is tangent to the worldsheet, \l i v k v = 0, i.e. when the unperturbed solution is invariant under the action of the 
relevant symmetry, but in such a case the ensuing strict conservation law is vacuous, since the current S{k} fJl will itself 



that can be seen from (37 
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be zero. In the case for which the unperturbed solution is invariant under the action of fc M one can however obtain 



a non trivial application of (40) by using the fact that in these circumstances if the first order dynamical equation 
( |33| ) is satisfied by a "real" displacement field £ M for a given 'real" gravitational perturbing field h pa , then the first 
order dynamical equation will also be satisfied by the Lie derivatives of these fields with respect to so these Lie 
derivatives will be utilisable used as the "virtual" fields in (fici|), i.e. we shall be able to take £ M = k v "V u ^ — (, y ^ v k v 
and h pa = k v V u h pa + 2h v ^ a \I p ^k v . 
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